Quantum critical phase diagram of bond alternating Ising model with Dzyaloshinskii-Moriya 

interaction: signature of ground state fidelity 
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We present the zero temperature phase diagram of the bond alternating Ising chain in the presence of 
Dzyaloshinskii-Moriya interaction. An abrupt change in ground state fidelity is a signature of quantum phase 
transition. We obtain the renormalization of fidelity in terms of quantum renormalization group without the 
need to know the ground state. We calculate the fidelity susceptibility and its scaling behavior close to quantum 
critical point (QCP) to find the critical exponent which governs the divergence of correlation length. The model 
consists of a long range antiferromagnetic order with nonzero staggered magnetization which is separated from 
a helical ordered phase at QCP. Our results state that the critical exponent is independent of the bond alternation 
parameter (A) while the maximum attainable helical order depends on A. 



I. INTRODUCTION 

One of the key features of strongly correlated electron sys- 
tems is their zero temperature behavior where quantum fluc- 
tuations have the dominant role to classify different novel 
phases which are separated by quantum phase transitions^. In 
this case, ground state (GS) is the sole candidate which re- 
ceives drastic changes at QCP driven by some external pa- 
rameters. Identification of QCPs is always a challenging 
task specially when quantum correlations diminish a single- 
particle picture. Within last couple of years different mea- 
sures of quantum entanglement have been proposed to be a 
new toolkit to detect and categorize QCPs^. Recently, ground 
state fidelity-the overlap of ground state at two slightly dif- 
ferent values of coupling constants-has attracted intensive at- 
tention as a proper quantity to signal QCP without the need to 
know the structure of phases close to phase transition 3,4 . 

An abrupt drop of fidelity in the vicinity of QCP is a con- 
sequence of an essential change in the structure of GS which 
is usually accompanied with a divergence of fidelity suscep- 
tibility. Meanwhile, it is rather difficult to obtain the GS of a 
strongly correlated many body system. However, recent suc- 
cessful attempts to calculate some measures of entanglement 5 
by quantum renormalization group (QRG) approach^ lead to 
a unified formalism^ to obtain fidelity in the thermodynamic 
limit (N —> oo) without knowing the GS of a large system. 
It suggests to implement QRG for obtaining the renormal- 
ization of fidelity for the bond alternating Ising model with 
Dzyaloshinskii-Moriya (DM) interaction. DM interaction 8,9 
which roots to the spin-orbit coupling is the antisymmetric 
super-exchange which leads to helical magnetic structures as 
the most likely candidates to host ferroelectricity 10 . 

In this article we study the quantum critical behavior of the 
one dimensional bond alternating Ising model with DM in- 
teraction in terms of renormalization of fidelity. The model 
represents an antiferromagnetic (Neel) long range order for 
small DM interaction while it undergoes via a continuous 
phase transition to a helical ordered phase. The ground state 
in the Neel phase is a product state which has essentially zero 
quantum entanglement and nonzero staggered magnetization 
while the helical phase poses a correlated quantum ground 



state with zero staggered magnetization. Within a classical 
picture the helical phase could be assumed as slightly rotating 
spins along the direction of chain. In terms of Landau the- 
ory of critical phenomena the staggered magnetization could 
be chosen as the proper order parameter. We have shown that 
the divergence in fidelity susceptibility (FS) at the QCP is an 
appropriate signature to find QCP which is more pronounced 
than the second derivative of ground state energyii. The pre- 
sented scheme makes us to find the critical points and its cor- 
responding exponents more accurately. 



II. QUANTUM RENORMALIZATION GROUP 

The Hamiltonian of bond alternating Ising chain with DM 
interaction on a periodic chain of N sites is defined 
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where Si is the spin- 1/2 operator at site i, A describes the rel- 
ative strength of the alternating coupling and J > shows 
the nearest-neighboring antiferromagnetic coupling. D is the 
vector of DM interaction which is considered in z-direction, 
i.e. D = Dz. To apply QRG&, the spin chain is decomposed 
to three-sites blocks (see FigJT} where the intra-block Hamil- 
tonian is H B and the inter-block one is H BB and their sum 
defines the whole Hamiltonian, H = H B +H BB (see the bot- 
tom part of FigJTJ . In this respect we have H B = Y^,f=i hf, 
where 
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FIG. 1. Decomposition of the lattice to 3-sites blocks. The intra- 
block Hamiltonian (H B ) and inter-block one (H BB ) are represented 
schematically in the bottom part. 



The block Hamiltonian (hf) is diagonalized exactly and the 
two lowest energy eigenvectors (l*^)) are kept to construct 
the embedding operator (Pj) to the renormalized Hilbert 
space, Pi = |*|) (fr | + |*7) (4| |. Here, | ff), | JJ.) represent 



the renamed states 



1*1) 



and I* j ), respectively in the renor- 



I*/) 



malized space to be considered as the new base kets. 
and 1*7) nave tne following expression (for odd blocks) in 
the original spin Hilbert space where | f) and | \.) represent 
the eigenvectors of S z operator at each site, 
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and Eq is the ground state energy of the block. The presence 
of bond alternation imposes to consider two types of blocks as 
depicted in FigQ] namely even and odd which are the mirror 
image of each other. For even blocks we find similar eigen- 
states by replacing a —c, b —> b and c — > —a. 

The global embedding operator is the direct product of the 
embedding operator of each block, P — ®^ Z Pi which gives 
the renormalized Hamiltonian by H' ren = P^HP 6 . The 
renormalized Hamiltonian is similar to the original one, EqQ] 
while the coupling constants are replaced with the renormal- 
ized one (denoted with ' ) as given in the following equations, 



J' = XJ, 



A' = A, 



D' 



Aab 2 c 



X 



-D, 



(6) 



where X = (1 - 2a 2 ) (1 - 2c 2 ). These relations, EqUJ de- 
fine the QRG-flow of our model which will be used in next 
sections and their features will be discussed in SeclIITl 



within the QRG iterations which hints to study the quantum 
phase transition by variation of D. The fidelity (f) associated 
to the ground states | *(£))) for a system of size N is defined 
by 



f(D,6;N) = (*(D-)\*(D+)), 



(7) 



where D± = D ± 5/2 and 6 is a small deviation around 
D. According to the renormalization group approach |*) = 
P|*W) in which P is the global embedding operator and 
l*^)) is the ground state of the renormalized Hamilto- 
nian. Thus, fidelity can be written in terms of renormalized 
ground state, / = (*( 1) (L»_)|P t (L»_)P( J D + )|*( 1) ( J D+)). 
A straightforward calculation shows that 

P}(D_)Pi(D+) = Ro(D_,D+)l, 
R a = [a(D_)a{D + ) + b(D_)b(D+) + c(D_)c(D + )},(8) 

for both even and odd types of blocks where X is the iden- 
tity operator. Therefore, the first QRG iteration leads to 

N 

f = R 3 (-D_,D+) x (*( 1 )( J D_)|*( 1 )(L> + )), where fidelity 
of the original model is expressed in terms of fidelity of the 

— 

renormalized one, i.e. / = i? 3 J' 1 '. It defines the renormal- 
ization of fidelity in terms of QRG. The QRG procedure is it- 
erated n-times to reach the renormalized system of N = 3™ +1 
and the ground state fidelity is expressed by 



f = d[RT T )(^(D.)\^(D + )), 



(9) 



i=0 



where Ri has the same expression as given in Eqj8]for Rq in 
which a, b and c are calculated at the i-th QRG iteration and 
(*(«)(D_)|*( n )(D + )) is the fidelity of a single block with 
three sites and n-times renormalized couplings. 

We have plotted fidelity (Eq|9]l versus D in Fig[2j-(left) for 
different chain length (N), 6 = 0.01 and at A = 0.5. By 
definition, fidelity is bounded by < / < 1 and an abrupt 
drop is a signature of quantum phase transition. We observe a 
sharp drop in Fig[2j-(left) for 0.8 < D < 0.9 which manifests 
that the ground state has encountered an essential change. The 
deep of drop is enhanced as the system size is increased which 
justifies an unfailing drop in the thermodynamic limit. This 
signature of quantum phase transition is more pronounced in 
the fidelity susceptibility (FS) which is the leading nonzero 
term in the expansion of fidelity and shows the change rate of 
fidelity, i.e. / = 1 - ^-FS + 0(5 3 ). Thus, FS is obtained by 
the following relation 



FS 



lim 2 

<5->0 



1-/ 

8 2 ' 



(10) 



A. Fidelity signature 

We implement the formalism introduced ii>£ to calculate 
the ground state fidelity in terms of quantum renormalization 
group. According to the QRG-flow (EqO, A does not run 



FigEKright) presents FS versus D for various system sizes, 
5 = 0.01 and at A = 0.5. A maximum appears in D = 
Dmax(N) which is increased by the size of system represent- 
ing a divergence in the thermodynamic limit. The position 
of D max (N) is exactly at the point where fidelity receives a 
drop. 
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FIG. 2. (Left) Ground state fidelity (f) versus Dzyaloshinskii-Moriya interaction strength (D) for different chain size and at A = 0.5. (Right) 
Fidelity susceptibility (FS) versus D at A = 0.5 for various system sizes N = 3 , . . . , 3 9 . Inset: the scaling behavior shows how the maximum 
of FS (Dmax) approaches the critical point (D c ) by increasing size. In both figures S = 0.01 which is the difference between two parameters 
(D+ — D- — 5) to calculate fidelity in the left figure and according to Eq.lllOt for the right figure. 



B. Scaling analysis 

It has been shown* that the fidelity susceptibility at the 
quantum critical point obeys a scaling relation. The scaling 
analysis for finite system size {N) states 

\D C -D max \~ N-i, (11) 

where D c is the quantum critical point, D rnax is the position 
of maximum in FS and v is the critical exponent governing the 
divergence of correlation length. The analysis of data of Fig|2} 
(right) is presented as an inset to this figure. It clearly verifies 
that the scaling relation EqQT|is satisfied with D c ~ 0.866 
and v ~ 2.17 for A = 0.5. Moreover, we have obtained 
similar behavior for different values of A (not presented here) 
where the critical point is found to be a function of bond alter- 
nating parameter, D c ~ y/l — A 2 while v ~ 2.17 is the same 
for all values of A. In contrast to what stated in Ref^ the crit- 
ical exponent which we got does not depend on A. It means 
that the whole phase boundary for < A < 1 belongs to the 
same universality class of a second order phase transition. 

III. SUMMARY AND DISCUSSIONS 

We have implemented the quantum renormalization group 
approach to study the zero temperature phase diagram of bond 
alternating Ising chain in the presence of DM interaction. To 
get a self similar Hamiltonian two types of blocks with 3-sites 
have been considered which leads to the QRG-flow equations 
of Eq|6] The QRG-flow tells that A does not vary within QRG 
procedure while D is renormalized. To get the phase diagram 
of the model we have calculated the renormalization of ground 
state fidelity which has been developed recently^. 



Fidelity as a geometric quantity— shows how much the 
ground state encounters an essential change by slightly mov- 
ing in the parameter space. Therefore, a sharp drop of fidelity 
versus a control parameter is a signature of quantum phase 
transition. The renormalization of fidelity obtained in Eq|9] 
in addition to QRG-flow, Eq|6] give the fidelity of our model 
for very large system sizes without the need to get the ground 
state. A clear drop of fidelity versus D in Fig[2j-(left) and 
consequently a maximum in the fidelity susceptibility, Fig|2j- 
(right), verifies the existence of a quantum phase transition at 
D c . We have applied the finite size scaling on the susceptibil- 
ity data presented in the inset of Fig[2]-(right) for A = 0.5 and 
generally got the critical phase boundary D c ~ \/l — A 2 . The 
phase boundary which separates the antiferromagnetic (AF) 
Neel order from a helical order is shown by the blue line in 
Figj3}(right). However, our scaling analysis gives a single 
value for the correlation length exponent v ~ 2.17 (indepen- 
dent of A) for the whole phase boundary which is in contrast 
\x£L. Although the presence of bond alternation breaks the 
translation invariance of the Hamiltonian it does not change 
the symmetry of the ground state which has already been 
spontaneously broken due to antiferromagnetic long range or- 
der. A comparison of our results with^ concludes that the 
bond alternation does not change the universality class of the 
model as far as A 7^ 1. 

We calculate the staggered magnetization (SM) and helical 
order parameter which is presented in Fig[3]-(left). Staggered 
magnetization is defined by 

1 N 

sM = -j2(-iywm, d2) 

i=l 

which can be expressed in terms of the renormalized ground 
state by replacing \^>) = P]^^). We use the projection of 
spin operators into the renormalized Hilbert space which fi- 
nally leads to 

SM = —SM {1 \ (13) 

o 
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FIG. 3. (Left) Staggered magnetization (SM, solid lines) and helical order parameter (Ch, dashed lines) versus D for A = 0.0,0.5,0.8. 
(Right) Zero temperature phase diagram in A — D plane. The blue curve represent the phase boundary between antiferromagnetic and helical 
ordered phases. The filled black circle denotes the first order tri-critical point and the red line shows a ferro-antiferromagnetic phase. 



where SM W is the staggered magnetization of the renormal- 
ized chain. A large number of iterations of EqQj] give the 
staggered magnetization in the thermodynamic limit. Simi- 
larly, the helical order parameter (Ch) is defined by the fol- 
lowing relation^ 

1 N 

i=l 

which can be calculated using the QRG-fiow. In this respect, 
Ch is expressed in terms of the helical order parameter in the 
renormalized model (C h ) by the following relation 

r _ MO) , r{0) Mi) 
C<°> = *£l£ , r(°)=-4a6 2 c. (15) 

o 

The above equation is iterated by replacing the couplings 
with the renormalized ones to reach the stable thermodynamic 
limit. We have plotted both staggered magnetization and he- 
lical order parameter versus D in Figj3j(left) for three values 
of A = 0.0,0.5,0.8 and N -> oo. For D = 0, SM is at 
its saturated value 0.5 and Ch = 0. The onset of nonzero D 
induces a helical order on the spins in the presence of anti- 
ferromagnetic order. The increase of D reduces the antifer- 
romagnetic order and enhances Ch (helical order). Exactly at 
the quantum critical point D c , SAI vanishes and remains zero 
for D > D c while Ch saturates to a finite amount which is 
less than its maximum attainable value. Similar behavior has 
been observed for all < A < 1 while the saturated value of 
Ch is slightly decreased by increase of A. 

To complete the phase diagram let us concentrate on the 
D = axis (Figj3]-(right)). At A = the model is sim- 
ply an antiferromagnetic Ising chain with Neel order | tl, tl 



, ■ ■ • , tl) and the ground state energy is Eq = —NJ/4. For 

< A < 1, there are two types of couplings (1 — A) J and 
(1 + A) J which are positive and induce the previous Neel 
order and ground state energy. Exactly at A = 1, the weak 
interaction (1 — A) J becomes zero and the spin model decom- 
poses to N/2 independent pairs of antiferromagnetically cou- 
pled spins. The ground state energy is still Eq = —NJ/4 
while the ground state is exponentially degenerate, namely 
2 JV / 2 . The degeneracy arises from the pairs which are de- 
coupled. For instance, | tl, t|, ■ ■ ■ , tl)> I It, tl, • ■ ■ , tl) 
and | |t, It, ' ' ' 7 It) are examples of different configura- 
tions which is possible for the ground state. The entropy 
(S) is proportional to ln(#of available states) which leads to 
S ~ ln(2 Ar / 2 ) = ^tt\ti2. This high amount of entropy at 
A = 1 is a signature of a phase transition which is actually of 
first order. Meanwhile, for A > 1 one of the interactions be- 
comes ferromagnetic, (1 — A) J < and the other (1 + A) J is 
still antiferromagnetic which totally leads to Eq = —NXJ/4. 
Thus, the derivative of Eq with respect to A receives a dis- 
continuity at A = 1. The ground state for A > 1 is either 
| U, It, tl, It, • • • > or | |t, U, It, U, " ■ ■)■ The first order 
tri-critical point is represented by the filled black circle in 
FigE}(right) and the ferro-antiferromagnetic ground state is 
denoted by the red line for A > 1 . 
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